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We obtain an expression for the energy of the density wave propagating in a multicom- 
ponent gravitating medium in the form well known from electrodynamics. Using the above, 
the possibility of "triple production" of the quasiparticles, or waves, with their energies 
summing up to zero, in a nonequilibrium medium is demonstrated. That kind of resonance 
wave interaction is shown to result in the development of an explosion instability. By the 
method developed in plasma physics, the characteristic time of the instability is evaluated. 

1 Introduction 

One would judge dark matter to be a natural object for the collisionless hydrodynamics 
of multicomponent gravitating systems, in which the particle movement is governed by 
the self-consistent gravitational potential. Up till now, its particles have been manifesting 
themselves solely in gravitational interaction, hence the lack of the evidence on the ques- 
tion of their nature (see for instance pQ). In spite of certain well-known achievements of 
the last years (like a bona fide Hubble photograph of 2004, of the two clasters colliding, 
which is strongly believed to have captured a ring of dark matter), the question of the very 
existence of the dark matter, raised by Fritz Zwicky in 1937 (see [2]), is still unaswered 

m. 

Development of the instabilities in a multicomponent system is determined by the 
hydrodynamical (or kinetic) parameters of the media composing it. Traces of structures 
or turbulences owing their appearance to this or that instability would be impossible 
to observe in the dark matter, but they should be present, say, in the gaseous disk. 
A list of instabilities capable of arising in the system "gaseous disk — dark matter" 
might come handy when one sets to drawing some conclusions concerning hydrodynamical 
(or kinetic) properties of the dark matter as a medium. Linear hydrodynamical beam- 
jeans instabilities of a multicomponent gravitating medium are described, for instance, in 
[I], [5J. If the media are sufficiently hot, the kinetic beam-jeans instability can develop 
instead of the hydrodynamical one. If the system is stable in the linear approximation (e. 
g., the large-scale instabilities are stabilized by rotation [6J), non-linear scenarios of the 
development of instabilities enter the game. 
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2 Energy of the density wave in a gravitating medium 



In this section, our aim is to obtain an expression for the energy of the density wave in a 
gravitating medium in the form convenient for further investigations. 

In [10], while deriving the formula expressing the energy of the density wave in an 
ordinary gas (i. e., self-gravitation amounting to effective zero), total energy of the unit 

2 

volume of a homogeneous medium is given as follows: E = pe + p being the density 
of the medium, e — the mass density of the internal density of the medium, v — its total 
velocity. Taking self-gravitation into account, one should add a term corresponding to the 
potential energy |/«/>, where ip is the gravitational potential. Apart from that following 
the standard procedure, for the density wave in the motionless medium we get: 

1 c 2 _ v 2 1 ~ 
E wave = -—p 2 + po— + -pip (1) 

Here the "tilde" symbol marks the perturbed values while the index "0" stands for the 
equilibrium, c s is the speed of the sound in the gas (we consider the isotermic case here, 
so that c 2 = p/p, p being the gas pressure). In a frame of reference where the medium 
has the velocity Vq, in the expression for the density of the energy of the wave one more 
term would appear, amounting to pv v. 

In the more general case of multicomponent gravitating medium with the plane-parallel 
relative motion of the components, the energy of the density wave per unit volume acquires 
the form: 

E wave = + PiVoiVi + + \pi^ (2) 

where the index "i" denotes the (both perturbed and pertaining to the equilibrium state) 
parameters of the component of that same number. 

The equations of gravihydro dynamics for the multicomponent gravitating medium 

have the following form: 

" & + divipiVi) = 



dt 

^ + (v,V)v 1 = (3) 



"i" being any integer from 1 to n, n — the number of the components of the system 
under consideration. 

Linearizing the equations and expanding them into Fourier series, for any given har- 
monic oc e - luJt + lkr we g e t 

_ p 0i k 2 ^ 

(w-kvoi) ~k 2 c 2 si 

kip (uj - kv oi ) , . 

v i = 2 \P) 

(w-kvoi) ~k 2 c 2 si 
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Substituting the above into the expression (j2J) for the energy of the density wave in a 
multicomponent medium, we get: 



E, 



p 0i u (uj - kv oi ) 



(6) 



wave 



((^-kv 0i ) 2 -^4) 2 



This is almost what we have aimed at. 

Now, the dispersion equation derived from the linearized equations of gravihydrody- 
namics: 



ujji being the jeans frequency of the component number i of the medium, defined in 
the usual way by uo 2 ^ = AnGpoi, G being the newtonian gravitational constant, and the 
physical interpretation of the value Eq will be discussed a few passages below. Direct 
calculations show that for the values to and k satisfying the dispersion equation (EJ), 



where g = — VV> is the perturbed gravitational field strength (for a given Fourier harmonic 



The structure of the expression ([8]) is fully analogous to the well-known results from 
classical electrodynamics [?] for the energy of the electromagnetic wave or the wave of 
the electric charge density in a dielectric medium. The difference rests in the sign (and in 
the presence, in the denominator, of the newtonian constant G), reflecting the difference 
in the nature of coulombian and newtonian interaction). 

The value Eq, analogous to the dielectric permettivity of a medium to the electric 
charge density waves, characterizes the responce of the density distribution in the various 
components of the gravitating medium to the perturbation of the gravitational potential, 
as it can be seen, for instance, from the expression for the perturbed density (j3J). 

So we see that the sign of the energy of a density wave in a multicomponent gravitating 

medium is determined by the sign of 

3 On the Possibility of Simultaneous Existence of the 
Waves of Positive and Negative Energy in a Mul- 
ticomponent Medium 

Consider the simplest case of a multicomponent medium — a case of a medium consisting 
of two streams of matter. Let both media composing the system have the same hydro- 
dynamical parameters (the isothermic sound speed c s and jeans frequency u 2 = 4irGp) 






(8) 



g 



ikip). 
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and move through each other with the velocities vq and — vq. For "one-dimensional" wave 
(k||v ) in such a medium, ([7]) takes the form: 



e = 1 + 1 + 1 = ( 9 ) 

(a; — kv ) — k 2 c 2 s {uo + kv ) — /c 2 c^ 

A sketch of behaviour of the left side of the dispersion equation (as dependent on 
frequency) is shown on Fig. 1. 
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Figure 1: The left side of the dispersion equation against the wave frequency 



If the curves lying in the lower semiplane intersect the straight line y — — 1, correspond- 
ing to Sq = 0, the dispersion equation has four real solutions (i. e., such a two-streamed 
system is stable in the linear approximation), two of them corresponding to the positive, 
and two — to the negative energy of the wave. (All the above holds true in the "super- 
sonic" case, when v > c s — in the other case the sketch on Fig. 1 would have looked 
differently.) 

It is easy to show that there exists a range of parameters in which such a system would 
be stable in the linear approximation. In this case, a wave propagating in the system can 
have either positive or negative sign in terms of energy. The dispersion equation Q can 
be reduced to biquadratic. Introduce a nondismensional constant r\ in such a way that: 

v 2 = (l + r])c 2 s 

1 The case of k_Lvo, as can be easily shown, corresponds to an ordinary jeans instability. But the waves 
of large wavelength {k 2 c 2 s < 2a->J) are excluded from out consideration anyway — by the limits of validity 
of our plane approximation as well as by the properties of real astronomical self-gravitating objects and 
their rotational geometry, in which the large-scale instabilities tend to be suppressed, see, for instance, 
i 
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The direct calculations show that the solutions u 2 of the dispersion equation ([9]) are real 
when 77 < 0. If 77 > 0, the solutions uo 2 are real whenever 

s — 2y/r)+l 

and amount to 



(uj 2 ) h2 = k 2 (2 + V )c 2 s - u) ± ^4(77 + \)Wc\ - 4{ V + l^c 2 ^ 2 + (10) 

Note that in the particular case k 2 c 2 = u 2 the solutions of the dispersion equation 
take the form: 



% = (2 + t^J 

,2\ _ m , ,2 



(t0 2 ) 2 = J]U)) 

We see that, in this case, if 77 > 0, the dispersion equation has four real solutions, two of 
them corresponding to the positive values of the wave energy and two — to the negative 
ones. 

Being interested in the solutions that are stable in the linear approximation (and thus 
corresponding to the ordinary (longtitudinal) density waves), we should now find k for 
which uj 2 is nonnegative. If 77 > 0, then for k satisfying the condition 



2 2 ^ 2 y/v + i + y/v 



all the four solutions to are real, two of them corresponding to the waves of the positive 
energy, and two — to the negative energy waves. 

Checking the case 77 = (^o = c s ), we see that the system is stable in the linear 
approximation if k 2 c 2 > uj 2 . If 77 < 0, the system demonstrates linear stabilities towards 
the perturbations for which k 2 c 2 > j^-. (One should remember that, when r] < 0, \r)\ < 1 
by definition.) 

4 The Resonance Three, or Non-linear Interaction 
Resulting in Explosion Instability 

Suppose that, being interested into all sorts of possible perturbations of our two-streamed 
medium, we stick to the range of linear stability. 

We now consider plane waves as quasi-particles. All further calculations will be carried 
out in the random phase approximation: we consider the waves (interpreted as quasipar- 
ticles) to be incoherent. The quantities necessary for evaluation of the probability of 
the event (resonance effect or any collision) are provided by the second correction to the 
gravihydrodynamical equations; the dispersion equation characterizing the dependence of 
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the frequency upon the wave vector we have already obtained from the first correction, 
or linear approximation. 

Quasi-particle collisions proceed with the conservation of energy and momentum. The 
important question is, can those conditions be held for a process sketched down on the 
scheme (Fig. 2) — triple quasi-particle production? 




Fig. 2: Triple Quasi-particle Production 

The conservation laws require that 

h + k 2 + k 2 = [ } 

As we shall see it below, such processes are indeed possible, and can be found even in 
one-dimensional model. Turn again to the waves with k||v . As it can be seen from the 
plot on Fig. 1, for each k, to the negative energy wave corresponds the pair of uj having 
lesser absolute value. Note by uj + the solution of the dispersion equation corresponding 
to the positive energy of the wave, and by uj~ - to the negative. We are interested in 
the "resonance three", meaning simultaneous production of, say, one "positive" wave and 
two "negative" ones. 



K) = fc?c*(2 + i7)-a;?- J4( V + l)kicj - A( V + l)k 2 c 2 s uj 2 + 



"1 



K) = k 2 c 2 (2 + rj)- ^-^4(77 + 1)^ -4(77 + 1)^0;] + ^ 

(O 2 = (A: i-A; 2 ) 2 cK2 + r ? )-^+ 

+ ^4(77 + l)(fci - k 2 ) 4 c A s - 4(77 + l)(h - k 2 yc 2 s cu] + 



(13) 



The first negative wave can be represented by a>j~, —k±, the second wave of negative energy 
- by uj 2 , k 2 , and the wave of the positive energy, to complete our "resonance three", will 
be represented by uo + = —uo^ — and the corresponding wave number k\ — k 2 . (The 
dispersion law that we have here is symmetric function with respect to uj (when k is 
fixed) as well as with respect to k with fixed uj.) In order to prove that such a "resonance 
three" indeed does exist, consider the case k 2 c 2 s » uj 2 . In the second approximation with 
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respect to ^f- we get: 

h = -y/rj + lk , LO X = W + 1 (VV + 1 - 1) k °s - 2kc a Vv+l 

k 2 = k , ^ = (^^+1-1)^-^^ (14) 
**=(V5FFT-l)* , ^3 = -^c s + J^pL 

Solving the equation uj\ + 0J2 + 003 = in terms of 77 gives the value 77 = | . For the short 
waves in the second approximation with respect to a small parameter we have presented 
"in quantities" the resonance three. The relative velocity of the two streams for the case 
under discussion must be 2v = 3c s . Thus the triple quasi-particle production, which can 
take place in a two-streamed gravitating medium, is described by the conservation laws 

m n 

Note by N k the number of quasi-particles carrying the 4-momentuno k (corresponding 
to the perturbations with the wave vector k and the frequency u>k, k and o;k being related 
to each other by the dispersion equation). Supposing N k » 1 one can write: 



O 7V T 

-ap = E Q{ki, k 2 ) [(N kl + l)(N k2 + l)(N^ k2 + 1) - N kl N k2 N^ k2 ] w 

(15) 



Q( k u k 2 ) [N kl N k2 + N k2 N_ kl _ k2 + N kl N. 



k 2 



Here Q{k\,k2) is the possibility of the triple quasi-particle production (or annihilation) 
of the quasi-particles with the momenta {k\, k 2 , — k\ — fe)- As it is noted in [9j, after 
the act of triple quasi-particle production is repeated quite many times, the number of 
quasi-particles of the three sorts that participate in the act can be considered roughly the 
same. Hence dismissing the indexes, we get: 

dN ,r2 

— = wN 2 
dt 

which leads us directly to the formula of the explosion instability: 

N(t) = (16) 

(Here iV is the number of quasi-particles of any of the three resonant sorts dependent 
on time, N - - the initial number of the quasi-particles. The coefficient w is to be 
interpreted as the possibility of the process, the value, inversely proportional to the time 
of development of the instability. In the next section we shall calculate it for the one- 
dimensional case.) 



2 Here we have put h = 1. 
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5 Evaluation of the Time of Development of the Ex- 
plosion Instability in a Two-Streamed Gravitating 
Medium 

The time to — {Now) -1 of development of the explosion instability (fl6|) is inversely pro- 
portional to the value w — the possibility of the resonance process responsible for this 
instability. (Naturally, as one can see from (jlBi) . w is the possibility per N, the number 
of quasi-particles, which is itself a dimensional quantity.) 

We shall use the method developed by certain authors for the decay instability in 
the random phases approximation in the beginning of 1960- ies. (See, for instance, [9].) 
Since the gaseous galactic disk, with its particles moving through the hypothetic halo of 
the dark matter, is in the turbulent state, the random phase approximation here seems 
adequate. 

The linear dispersion equation was obtained in the first approximation with respect to 
the amplitude; considering the case of weak nonlinearity, one should turn to the second 
approximation. 

Consider the Fourier component of the Poisson equation correcponding to the momen- 
tum k and frequency uj. For the sake of convenience introduce the 4-dimensional wave 
vector k = (w,k). Whenever the inverse is not clearly stated, speaking of the 4- vector k 
we would always have in mind k = (c^k, k), where = uj(k) is the solution of the disper- 
sion equation £o(u>, k) = 0. By the above relation a manifold in the 4-dimensional space is 
defined; the 4-dimensional integration appearing in the further calculations is performed 
over this manifold. Below, among the indexes of the Fourier harmonics, 4- vector k as well 
as the normal 3-vector k, will be used. 

- k%(^, k)^ fc = 4vrG( P ; 2 fc ) + pg») (17) 

Here, as well as in the previous sections, Eq is the gravitational analogue of the dielectric 
permettivity. Index (2) over the symbol of density p indicates the correction of the 
second approximation. In the linear case Eq(u, k) = would have held; if we present the 
small nonlinear correction to the frequency as Sujn + ijN, then, as is easily seen, for our 
approximation 

(buj N + I'Jn) 

Multiply both parts of the equality by some k'-component of the perturbed gravitational 
potential ?/v and then average it over phases: 



*o(o;,k)« — 



-{buj N + i'y N ) — 



k 2 I k 5(k + k') = AnG ( (pg + pS) V* 



Here Ik is the spectral intensity (< ipk4>k' >— Ik${k + k') in the random phases approxi- 
mation). 
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The energy of the plane density wave is 



W k = - J- ^1 

k 8nG duo 



\9k\ 



|2 



Since 5^ = — (VV>)fc — — ^k^fc an d £o(^k) = 0, for the number of quasi-particle with the 
4-momentum k one can write: 



N = Wk = l_ deo 

uJk 8ttG duo 



k 2 4 



In the expression for the nonlinear frequency shift Su^f+i'y^ the real part can be neglected 
compared to — it corresponds to the renormalization of the linear dispersion law. 
The imaginary part 27^ corresponds to the nonlinear spectre evolution. Note that for 
ip ~ e -*(^+*7Jv) number of quasi-particles will be ~ e 27jv *, so that multiplification by 
7jv for the function Nk(t) is equivalent to the action of the operator of the differentiation 
over the time: 7jv ~ §J^- Taking all this into account, note that the equation (fT7j) 
transforms into: 

dN, 



dt 

On the other hand, the structure of the time evolution of the number of quasi-particles 
is determined by the equation (fl5|) . or, in the integral form, 



dN k 

dt 



J io(k, q)<5(^ k + cu q + c^_ k _ q ) (jV q M k _ q + iVk-^Vq + iV k iV- k - q ) dq (19) 



To make the long calculations short, we now set to follow the evolution of only one 
term of the sum under the symbol of integration, say, iV q iV_ k _ q . We shall sort of exstract 
it from the second correction to the density; the analogous expressions for the rest two 
terms will be easy to construct; besides, the coefficient (i. e. the possibility in which we 
are interested) preceding them in the expressions like [19] is the same. We take the second 
correction to the density in the following form: 



Pak = J < q (ipql/l-k-q ~ (^-k-q)) dq (20) 

Here a = 1, 2 is the number of the corresponding component of the medium, is the co- 
efficient to be calculated a few passages later, from the Euler equations. The 3-dimensional 
indexes k, q are used instead of the 4-dimensional ones because three components of the 
wave vector are sufficient to fully define the coefficient v. The average value (ipgip^k-q) 
corresponds to the evolution of the phone and should be subtracted. 
Substituting ([20]) into 02]), we get: 




^S(k + k') = -i I ( > ; < q 1 (W {%^-k-q- < M-k-i >)) dq (21) 
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Substitute now the Fourier-component ipk> by its nonlinear expansion (through the 
density from (ITT)) ): 



4 ] 



< Ipgilp-k'-q, >) dq' 



in the course of this substitution we would have to perform the averaging procedure. The 
randomness of the phases results in reducing of all the correlations to binary, so that: 

{ipq^-k'-q'lpqlp-k-q) =< ^q'^Pq >< Ip-k'-q'^-k-q > + < ^q'^-k'-q' >< ^qlp-k-q > + 
+ < VYV'-fc-g >< ^-k'-q'^q > 

The middle term here vanishes by the virtue of (I2TI) . and the rest two terms are symmet- 
rical with respect to replacement q' — > —k' — q', so that one of them can be doubled and 
the other one written off: 



(VyV'-fc'-g'V'gV'-fc-?} = 2J,I_ Jfe _ 5 % + q')5{k + k') 



2 ■ (8nG) 2 N q N_ k _ q 



(k + q) 



2 de 
du 



,2 ^£o 
du 



q 



-i 



X 



x 5(q + q')5(k + k') 
Using an approximate equality in terms of generalized functions: 



^K^k') 



\u — u 



de 











I du 


<"=<W 



-1 



illS(u' — u' k t) 



and substituting all that into (I2TT) . performing integration over <ig' and dk', for the evolu- 
tion of the number of quasi-particles (tracing the transformations of one of the first two 
terms of ( IT9l) — iV q A^_k-q — we end up with the following expression: 



dN k 
dt 



IT 



dq(8nG) 3 (J2 



< q ) 2 X 



k 2 ^£o 

du 

X N q N-k-q 



(k + q) 



de 
du 



q" 



du 



-i 



x 



:22i 



Here r is the dimensional coefficient, by the order of magnitude similar to L/vq] 
t = j^prg J S(u' — c^kO^k'. Thus we see that the expression for the probability has the 
following structure: 



w 



(87rG) : 



7TT 



1A 



k,q 



du 



k-q) 



eta 



q~ 



du 



(23) 
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Now one has to determine the coefficient z/£ q given by relation (I20I) . The second 
correction to the continuity equation is: 



^ + pWdi«v« + (vWV)pW + po^^vi 2 ) + (v 0a V) P L 2 ) = (24) 



The Euler equation in the second approximation: 



<9v 



(2) 



Of 



(vi 1) V)v« - (vWV)v. 



(0)VW( 2 ) 



(25) 



Fourier-expansion for (a?, k)-component of the continuity equation gives: 



"pgl = / dqpg ((k - q), vgi _ q ) + | dq (q, v«_ q ) fl& + Po. (kv« ) + pi 2) k (kv 0a ) 

(26) 

Rewriting the same relation in more symmetrical form: 



pi 2) k (^ - kv 0Q ) = - J dq [pW(k, < k _ q ) + P^- q (k, vgyj + p a(k, < k 

The values of the linear approximation are related in the following way: 

kv (1) 

(1) _ KV Q,k 

Pa,k - P«0" 



(i) 



,( 2 ) 



(27) 



- kv 0Q 

Fourier transform applied to the Euler equations gives: 



(- " kv te )v2 = J / ^ Vg (q, V«_ q ) + v£ q (k - q, V«) 



r (l) 



:28i 



In order to obtain the result in the analitical form, we now turn to 1-D case. As in 
previous sections, let p i = P02 = Poi ^oi — —^02 = v o- Below the index \\kn, marking 
the Fourier component, will mean now the one- dimensional "wave number", not the 4- 
momentum. 

From the Fourier transform of the linear approximation one easily obtains: 



(i) _ k(u k ~ kv ) 

Vl > k ~ (uj k - kv ) 2 - k?c\ 

(i) k 2 p 
Pik 



(cu k - kv ) 2 - k 2 c 2 s 

Changing vq — > — vq and (1) to (2) in the above expression, we get the similar relation for 
the second component. 

Substituting (I28I) into fl27|) and taking into account the above expressions for the 
Fourier components of the velocity and density in the linear approximations as related to 
the corresponding Fourier component of the amplitude of the gravitational potential, also 
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having in mind that v x _ kq = v\ q due to the symmetry of the initial conditions, we obtain 
the coefficient: 

„2 * kq{k + q) v7 ,2,3/ \ / 9Q n 

= 2"° x x H (29) 

Here (in the above) for each x = k, q, —k — q the following notations have been 
introduced: 

A i = (^) 2 - X 2 cl 

and 

7?' s (±) = ?n^_ ff n^ + (-* - ?)njn^ + 

For the first component of the medium (the stream moving with the positive velocity 
Vq) the corresponding coefficient v\ can be obtained from v\ by replacement v — > — 1> . 
In the above passages, performing the Fourier transform, we have not usually mentioned 
the dimensional coefficient; it is easy to see that the latter will be important only in 
calculation of the coefficient r, which in one-dimensional model is equal to: 

L being the characteristic scale of the problem, the linear size of the system in the direction 
parallel to the velocity of both streams. Calculating this expression for k' = k, Uk> = u>k, 
we get: 

£ ( A t) 2 V- + (^)^+ 

2tt (-lovo + k V c 2 s )(A+) 2 + (uvo + k V c 2 s )(A-) 2 

Thus the probability is, effectively, found. The expression for the probability even in 
one-dimensional case looks somewhat awkward: 

= ^rr (T**(-)A+At k _ q (tot k )* + T^(+)A q AZ k _ q (toZ k ) 2 ) 2 (^) 2 (A+) 2 
W 4p X (nt) 4 (^) 4 D k D q D_ k _ q [ } 

Here for all x — k, q, —k — q the following notation is introduced: 

The expression obtained above, in principle, provides a possibility to evaluate the time 
of development of the instability, if the hydrodynamical parameters of the problem (i. 
e. density, speed of sound) are known and the "resonance three" are determined. In 
our case, since the two streams going through each other in the opposite directions are 
hydrodynamically identical, and the "resonance three" has been obtained, all we have to 
know to evaluate the time of development of the explosion instability is the density and 
the speed of sound. 
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For the resonance three determined in one of the previous sections, if the characteristic 
size L is taken to be the jeans size of one of the components, then the value w, calculated 
according to (!30il . acquires the following structure: 



4ttuj 2 

w - 



p c 2 s (3 

(3 being the relation of the jeans wavelength to the wavelength of the perturbation in 
question. Supposing that in gravihydrodynamics, like in plasma physics, the degree of 
the turbulence can be evaluated by the parameter << 1, where W is the density 
of the noise energy, we get the evaluation for the time of development of the explosion 
instability: 

t ins ~ a(k)- 



noTujj 

a(k) being a dimensionless coefficient which is the less the greater the turbulence degree 
of the medium is. If W(k) has a pointed peak then for the value of k corresponding 
to that maximum a(k) << 1. Since, obviously, the second factor in the expression for 
t ins has the same order of magnitude as the inverse jeans frequency (say, 10 15 — 10 16 
s for the interstellar gas), that would mean that, for those k, the instability would be 
able to develop. It is important to have in mind that what we deal with is a nonlinear 
process, achieving the infinite amplitude in finite times (ideally speaking). Comparing 
linear exponential instabilities we take the definitely fastest one in terms of increments; 
to say that one of the instabilities is fast enough to develop against the background of 
the other, we require a strong inequality 71 >> 72 (71, 72 being the increments). But 
when one of the instabilities in question is explosion instability — the simple inequality 
is enough to conclude that it is fasteiH 

Returning to the real astronomical situation, it should be noted that the matter of 
the gaseous disk and the dark matter can hardly be considered hydronamically identical. 
To take the hydrodynamical difference between the two streams into account should not 
change dramatically the expression for the probability: each coefficient u^ q , a = 1, 2, will 
be expressed through the equilibrium density of the corresponding component, and the 
auxiliary constructions like would get additional indexes, so that: 

A+ = (u x + X Vo) 2 ~ k 2 c 2 sl 

A x = ( u x ~ XVo) 2 - k 2 c 2 2 

Yet the dispersion equation would not be analitically solvable any more; numerical calcu- 
lations would be required. 



3 In our case the actual characteristic sizes are not, in fact, jeans wavelengths: they are determined by 
the conditions limiting the applicability of our plane approximation. In the real cosmic space geometry, 
however, the large-scale instabilities are stabilized by rotation. 
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6 Conclusions 



In the present work, we have obtained the expression for the density wave in a gravitating 
medium, analogous to the corresponding formula in plasma physics. The analogy has 
been toyed with by the theoretists of both fields, but seems to have never been clearly 
formulated and demonstrated in the general case. The important role in this parallel 
between the two regions of physics is played by a physical value describing the responce 
of the distribution of the density of the gravitating matter to a change in the gravita- 
tional field. In electrostatics and electrodynamics the analogous value is the dielectric 
permettivity of the medium. 

Having thus legitimized certain analogies between the plasma physics and physics of 
gravitating media that have been noted before, we concentrated on development of one 
of those analogies, related to the possibility of development of the explosion instability in 
a two-streamed gravitating medium. 

We have shown that, just as in a two-streamed plasma, in the two-streamed gravi- 
tating medium simultaneous propagation of the waves of positive and negative energies 
is possible. On the example of two streams moving through each other in opposite di- 
rections and having the same hydrodynamic characteristics it was demonstrated that a 
multicomponent gravitating medium can provide the background for development of the 
explosion instability. 

The explosion instability is faster than the linear, exponential instabilities; the ampli- 
tude of the perturbation reaches infinity within finite time. Using the standard method 
developed in plasma physics, we have demonstrated the way to calculate this time for 
one-dimensional model. 

Stabilization of the explosion instability is achieved through restructuralization of the 
dispersion equation. The dependence of the frequency of the perturbation on the wave 
vector that we have used in our calculations was obtained from the linear approximation. 
For large enough amplitudes, this dependence changes its form, the "resonance three" 
do not resonate any more, so that triple production of these quasi-particles becomes 
impossible. 

Classification of various scenarios, composition of a librarian's catalogue of all pos- 
sible instabilities has been long recognised as an important task in plasma physics. In 
graviphysics, the instabilities responsible for the most part of the observed phenomena 
and structures in the Universe probably deserve no less attention on the part of physicists 
and astronomers. 

The authors would like to express their gratitude to A. L. Barabanov, K. V. Chukbar 
and A. M. Fridman for fruitful discussions. 
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